Asymmetric double-well potential for single-atom interferometry 
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We consider the evolution of a single-atom wavefunction in a time-dependent double-well inter- 
ferometer in the presence of a spatially asymmetric potential. We examine a case where a single 
trapping potential is split into an asymmetric double well and then recombined again. The inter- 
ferometer involves a measurement of the first excited state population as a sensitive measure of the 
asymmetric potential. Based on a two-mode approximation a Bloch vector model provides a simple 
and satisfactory description of the dynamical evolution. We discuss the roles of adiabaticity and 
asymmetry in the double-well interferometer. The Bloch model allows us to account for the effects 
of asymmetry on the excited state population throughout the interferometric process and to choose 
the appropriate splitting, holding and recombination periods in order to maximize the output signal. 
We also compare the outcomes of the Bloch vector model with the results of numerical simulations 
of the multi-state time-dependent Schrodinger equation. 

PACS numbers: 39.20.+q, 03.75.Dg, 03.75.Be 



I. INTRODUCTION 



The evolution of a quantum system in a double-well 
potential has been the subject of numerous theoretical 
studies. These include treatments of Josephson-like os- 
cillations Q, 0, dynamic splitting Q and interference 
of Bose- Einstein condensates (BECs). The interpre- 
tation of these effects is based on the approach in 
which interference patterns are seen to evolve as a re- 
sult of successive boson measurements which do not iden- 
tify the originating condensate. The production of cold 
atoms and BEC in microtraps on atom chips H, 0, 13 
and in micro-optical systems |9j has stimulated a great 
interest towards novel implement at ions of atom interfer- 
ometers [13, El HI 13 Q El that are based on the 
use of double- well potentials. Double- well atom interfer- 
ometers (DWAI) of both the single-atom and the BEC 
varieties are well suited to implementation on an atom 
chip. Here microfabricated structures allow a precise con- 
trol on a sub- micron scale over the splitting and merging 
processes. The key processes of splitting JTfJUJlllS an d 
merging fTfjj of cold atomic clouds and even interference 
of a BEC after splitting in a double well pfl El 15^ 
have been already demonstrated. Although the imple- 
mentation of a DWAI using a BEC can lead to a \/~N- 
fold enhancement in precision measurements j23l |. phase 
diffusion [24| associated with mean field effects is of con- 
cern |2fJ. Double- well interferometry with a single atom 
can allow a longer measurement time and in this regard 
has a potential advantage. An on-chip single-atom inter- 
ferometer can be integrated with the source of atoms in 
a ground state - the Bose-Einstein condensate - and be 
used for sensitive measurements of gravitational fields. 
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DWAI may also be applied to measure collisional phase 
shift induced by the atom-atom interaction, which is use- 
ful for quantum computation processes |25j. 

Two proposed schemes of a single-atom DWAI involve 
time-dependent transverse |lfj| and axial splittings 
of a trapping potential. An atom is initially prepared in 
the ground state of a single symmetric trapping poten- 
tial, which is then split into a symmetric double well. A 
spatially-asymmetric potential is then applied and non- 
adiabatic evolution leads to transitions between ground 
and excited states. The asymmetry is then switched off 
and the double well is recombined into the original po- 
tential. The population of the excited state measures 
the effect of the asymmetric potential. The DWAI can 
be considered as a quantum-state Mach-Zehnder inter- 
ferometer where the evolution of the quantum state via 
the two separated wells is analogous to the propagation 
of an optical field via two pathways. 

However, these treatments ignore the effect of asym- 
metry during the splitting and merging stages. In real- 
ity asymmetry is always present and could be the result 
of imperfect horizontal splitting (introducing a gravity- 
based asymmetry), external spatially- variable magnetic 
and electric fields or different left and right trap fre- 
quencies. We show that the presence of small asym- 
metries have dramatic consequences on the interferomet- 
ric process. We have produced a simple model in terms 
of a Bloch vector evolution that enables us to consider 
a splitting-holding-merging sequence involving a double 
well and takes into account the presence of asymmetry 
at all stages. The two key parameters are the energy gap 
between the lowest two states of the symmetric compo- 
nent of the trapping potential and an asymmetry param- 
eter, which is related to matrix elements of the asymmet- 
ric component of the potential. Non-adiabatic evolution 
only occurs during the splitting and recombining stages 
when the torque vector changes much more rapidly com- 
pared to the Larmor precession of the Bloch vector. It is 
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important that the torque vector remains constant during 
the holding stage, and that this period is long compared 
to the splitting and recombination times. In this case 
the final excited state population is a sinusoidal function 
of the holding time, with a period determined via the 
asymmetry parameter. 

In this paper we consider the dynamics of a single 
atom in an asymmetric DWAI, with the basic theoret- 
ical framework being outlined in Section II. Using a two- 
mode approximation we develop a Bloch vector model for 
time-dependent DWAI (Section III) that provides an ad- 
equate description of the dynamics of the splitting, hold- 
ing and recombination processes in the presence of an 
asymmetric component, and is then used (Section IV) in 
describing the interferometric process. The validity of the 
two-mode approach is explored in Section V through the 
comparison of predictions of the Bloch vector model with 
the results of direct numerical simulations of the time- 
dependent multi-mode Schrodinger equation. A discus- 
sion of results follows in Section VI and includes a novel 
scheme to measure the population of the excited state. 
Theoretical details are dealt with in the Appendix. 



II. THEORETICAL FRAME 

In general, the evolution of a single atom in an inter- 
ferometer must be described using a three dimensional 
quantum treatment. However, for a system of cylindri- 
cal symmetry (as is present in typical atom chip experi- 
ments) it is possible to ignore excitations associated with 
the two tightly confined dimensions, as long as the dy- 
namics throughout the process is restricted to the dimen- 
sion of weak confinement (longitudinal splitting) . In this 
system it is possible to reduce the quantum treatment to 
that of a one dimensional problem. 

We consider the one-dimensional evolution of a single 
atom system due to a time-dependent Hamiltonian H(t) 
that can be written as the sum of a symmetric Hamilto- 
nian Ho(t) and an asymmetric potential V as (x) 
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where a specific form for the symmetric potential Vq is 
chosen 26] . The Hamiltonian and other physical quan- 
tities have been written in dimensionless quantum har- 
monic oscillator units associated with atomic mass m and 
angular frequency loq. With the original quantities de- 
noted by primes we have 
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The dimensionless Hamiltonians, potentials and energies 
are obtained by dividing the original quantities by Tiujo- 
V as will be taken as a linear function of x. 

The symmetric potential depends on a time-dependent 
splitting parameter /3, whose change from zero to a large 
value and back to zero again conveniently describes the 
splitting, holding and recombination processes with peri- 
ods T s , Th and T r respectively. For zero /3 the symmetric 
potential involves a single quartic well. When it is large a 
double harmonic well appears with a separation between 
minima of 2y/2f]. For zero [3 and for large x the symmet- 
ric potential approximates that for a quantum harmonic 
oscillator with frequency too and mass m. Key results in 
the paper would still apply if other suitable forms for the 
symmetric potential are used. 

We denote the eigenvectors of H as \<fii) and their en- 
ergy eigenvalues as Ei, where i = 0, 1, 2, .. and Ei + i > Ei. 
The corresponding quantities for the symmetric compo- 
nent of the Hamiltonian, Ho, will be denoted \Si) and 
Esi- Both sets of eigenvectors are orthonormal, and 
all energies and eigenvectors are time dependent. Ho 
is symmetric and the ground state \So) is symmetric and 
denoted as IS"), whilst the first excited state is an- 
tisymmetric and denoted as \AS). Their energies are 
denoted E$ and Eas- The one-dimensional nature of 
the system allows real eigenfunctions <pi(x\ Si(x) to be 
chosen. In this case the geometric phase |27jj is zero. 

We can illustrate the general behavior of the lowest few 
energy eigenvalues (Eo, E\, E2, E3, ..) of the total Hamil- 
tonian H as the splitting parameter is increased from zero 
to a large value and back. At the beginning and the end 
of the process when [3—Q the energy eigenvalues are well 
separated. Here the effect of asymmetry V as is small and 
the eigenvalues and eigenvectors resemble those for the 
symmetric Hamiltonian. When the splitting parameter 
increases and the trapping potential changes to a double 
well, pairs of eigenvalues (Eo and E\, E2 and E3, etc) 
become very close. At this stage the quantum system 
is very sensitive to the presence of V as which breaks the 
symmetry, allows transitions between |^o) and \4>i) to 
occur and causes the eigenvectors |</>o) and |0i) (as well 
as |(/)2) and |</>3)) to be localized in the individual wells 
in this far split regime. 

Initially, the atom is prepared in the lowest energy 
eigenstate \cf>o) of the single well. Transitions to excited 
states are suppressed if the time scale for splitting and 
recombination is much longer than the inverse frequency 
gap between the relevant states. The energy gap be- 
tween Eo and E2 is always larger than the gap between 
Eo and E\, and by choosing appropriate time scales we 
can isolate the two lowest energy eigenstates (\<j>o) and 
\cj)\)) from higher excited states (|</>2), ^3), etc), but still 
allow for transitions between the two lowest energy eigen- 
states to occur. As a consequence the dynamics of the 
DWAI can be treated under the two-mode approxima- 
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tion, in which only the two lowest energy eigenstates of 
the total Hamiltonian H and the symmetric Hamiltonian 
Ho need to be considered. In this case the first excited 
state probability (a measurable quantity) can vary from 
zero to one. A proposal for measuring the excited state 
population is outlined in Section VI. 

Using the two-mode approximation expressions for the 
lowest two energy eigenvalues (Eq, E\) and eigenvectors 
i\4>o)i l^i)) f° r the Hamiltonian H will be obtained. A 
standard matrix mechanics approach will be used, but 
instead of using the symmetric potential energy eigen- 
vectors \S), \AS) as basis vectors, we use the left, right 
(L-R) basis vectors \L), \R), which are defined by 

\L) = -^(|S) + |AS», (4) 

\R) = ±(\ S )-\AS)). 

The states \L), \R) are orthonormal and for large j3 
correspond to an atom localized in the left or right well, 
respectively. However, even for a single well the L — R 
basis is still applicable. The matrix for the Hamiltonian 
H in the L — R basis is given by 

where the order of the columns and rows is L, R and we 
define the convenient real quantities 

A = E AS -E S , (6) 
V as = (R\vZ\R)-(L\vr s \L) 

= -(<S|t£|AS> + <AS|t£|£>). (7) 

The derivation of the Hamiltonian matrix uses the 
symmetry properties of \S), \AS) and the reality of 
the eigenfunctions. The total energy for the symmetric 
Hamiltonian is given by Es + Eas, and the energy gap 
is given by A . The quantity V as describes the asymme- 
try of the system, and would be zero if the Hamiltonian 
was symmetric. The second equation relates V as to off- 
diagonal elements of the asymmetric contribution to the 
Hamiltonian, indicating its role in causing transitions be- 
tween the eigenstates \S), \AS) of the symmetric Hamil- 
tonian. The Hamiltonian matrix JSJ is analogous to that 
for a two-level atom interacting with a monochromatic 
light field 28]. The symmetric Hamiltonian transition 
frequency Ao is analogous to the Rabi frequency, whilst 
the quantity V as is analogous to the detuning. 

The energy eigenvalues for the total Hamiltonian H 
are obtained from the determinental equation as the 
eigenvalues of the matrix [H ] L ~ R , and are given by 

E = \{E S + E AS )~ I A, (8) 
Ei = ^(E S + E AS ) + ^A, 



where the quantity A gives the energy gap for the total 
Hamiltonian H and is defined by 

A = yj Al + V* =E l - E . (9) 

In terms of the laser-driven two-level atom analogy, A 
would be analogous to the generalized Rabi frequency. 

The orthonormal energy eigenvectors for the total 
Hamiltonian H are given by 

I0o> = \L) + \J 1 -^- \R), (10) 

Ifc) = yjl^\L)-^^\R), 

where the effect of asymmetry is now represented by the 
variable 

V = -^. (11) 

A v ' 

On substituting for \L), \R) the eigenvectors for the to- 
tal Hamiltonian can be related to those for the symmetric 
Hamiltonian Hq . At the beginning and end of the inter- 
ferometer process we find that the asymmetry parameter 
V as is small compared to the symmetric energy gap Ao. 
For V <§; 1 the eigenvectors \4>o), \<j>i) become similar to 
IS") and \AS) respectively. For V as > A (V ~ 1), the 
eigenvectors \4>o), \(j>i) become equal to \L), \R) respec- 
tively, the localized eigenvectors for the separate wells. 

The behavior of the quantities A, V as and A as the 
splitting parameter (3 is changed is shown in Fig. ^ for 
the case where the asymmetric potential V as (x) varies lin- 
early with the coordinate x, specifically with V as =0.02x. 
The symmetric energy gap Ao becomes close to zero for 
j3 > 4 and then the actual energy gap A is approxi- 
mately given by V as . The energy eigenfunctions (j)o(x) 
and 4>i(x) for different splitting parameters (3 are de- 
picted in Fig. [3 again with V as =0.02a;. The behavior 
of the potential V(x) = Vq(x) + V as {x) is also shown. 
For small (3 [Fig. Ufa)] the potential is a single well and 
the eigenfunctions are approximately symmetric and an- 
tisymmetric. For larger f3 [Fig. |2{c)] the potential is a 
double well, which still appears to be symmetric. How- 
ever, even with a small asymmetry in the potential the 
eigenfunctions are no longer symmetric and antisymmet- 
ric, but instead are each localized in separate wells. This 
sensitivity of the eigenfunctions to a small asymmetry is 
critical to the performance of the present interferometer. 

III. BLOCH VECTOR MODEL 

We can express a general time-dependent normalized 
state vector |<3>) as a quantum superposition of the states 
\L) and \R) 

\*(t))=C L \L)+C R \R). (12) 
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FIG. 1: Energy difference A between ground and first excited 
states (solid line) for V as =0.02 function of the splitting 

parameter fl. Dotted line - energy difference Aq for symmetric 
Hamiltonian, dashed line - asymmetry quantity V as . 

Our interferometer will be described in terms of the 
Bloch vector and its dynamics determined via Bloch 
equations. We now introduce Pauli spin operators and 
the Bloch vector. Time-dependent Pauli spin operators 
Oa (a = x, y, z) are defined in the Schrodinger picture 

S£ = (\R) (L\ + \L) (R\), 

<T y = \ {\R)(L\ \L){R\), (13) 

a~ z = (\R) (R\ \L) (L\). 

From its matrix representation in the L — R basis 
(jSJl, the dimcnsionless Hamiltonian operator H in the 
Schrodinger picture can be expressed in terms of the 
Pauli spin operators as 

h = -(o, i + n x o^ + n y -^ + n z s^), (14) 

where 

fio = {Es+Eas), (15) 

Sl x = — Ao, = 0, Vl z — V as - 

It is convenient to introduce a so-called torque vector, 
defined as (l = (fi x , ft y , D, z ). 

The Bloch vector is defined to have components which 
are the expectation values of the Pauli spin operators 
Oa in the quantum state I'F). These components will be 
denoted as a a - Hence in the Schrodinger picture 

<r = I *(*)) {a = x,y,z). (16) 

The Bloch vector is defined as a = (a x , a y , a z ). The 
Bloch components are bilinear functions of the ampli- 
tudes Cl and Cr. 
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FIG. 2: Stationary eigenfunctions of the ground state 
(dashed-dotted line) and the first excited state (dotted line) 
for different values of the splitting parameter /3 = 1 (a), 2.5 
(b) and 5 (c). The potential V(x) is shown as the solid line. 



The evolution of the DWAI system is now described 
by a set of real variables a Xl a y , a z and each of these 
variables has a certain physical meaning. The component 
a z is a measure of the imbalance of the atomic population 
of the localized states \L), \R). The component a x is 
a measure of the atomic population imbalance between 
the delocalized states IS*), 1^45), as can be seen if the 
quantum state is expanded in the symmetric basis. For 
a x — +1 all the population is in the symmetric state IS"), 
for a x = —1 it is all in the antisymmetric state \AS). 
It is thus a measure of the excitation of the first excited 
state in the unsplit trap regime. 

Equations for the components of the Bloch vector can 
be obtained from Heisenberg equations for the Pauli spin 
operators. The derivation must take into account the 
present situation where the Pauli spin operators are ex- 
plicitly time dependent since the basis vectors |L), \R) 
change with time. This differs from the standard sit- 
uation of time independent basis vectors [2^, . How- 
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FIG. 3: Evolution of the Bloch vector a and the torque vector 
Q at different moments of the splitting stage: (a) - at the 
beginning when Ao 3> V as and H ~ (— Ao, 0, 0) (b) - when 
A = V as , (c) - when V as > A and Q « (0, 0, V as ). 

ever, the additional term in the Heisenberg equations can 
be shown to contribute zero to the Bloch equations due 
to the two eigenfunctions in the symmetric basis being 
real and having opposite symmetry (see Appendix). The 
Bloch equations are given by (dt = d/dt) 

dt&x = ~V as <J yi 

d t cr y = V as a x + A cr z , (17) 
d t a z = -A Q a y , 

and can be solved numerically using the Runge-Kutta 
algorithm. In vector notation the Bloch equations are 

d t S = Vlxa. (18) 

This form of the Bloch equations is a direct conse- 
quence of the equivalence of the two-mode double-well 
interferometer to a spin 1/2 system. The Bloch vector 
precesses at the Larmor frequency around the torque vec- 
tor, which in detail is 

d = (-A , 0, Vas). (19) 

If there is no asymmetry, the x component of the Bloch 
vector remains unchanged, whilst its component in the 
y — z plane just rotates about the x axis (Fig.[3J). 

IV. MODEL OF A SINGLE-ATOM 
DOUBLE- WELL ATOM INTERFEROMETER 

In the single-atom interferometer under consideration 
the atom is always located in a trapping potential, which 



changes from a single well to a double well - which in 
general is slightly asymmetric - and back again to the 
original single well. The interferometer is used to mea- 
sure the effects of this asymmetry, the cause of which may 
be of measurable interest (e.g., as in a gravity gradiomc- 
ter). The atom is initially in the ground state \<po (0)) of 
the original unsplit potential, and as V as is then small 
compared to Ao, \4>q(0)) is then approximately the same 
as |5(0)). The population of the excited state at the end 
of the recombination process is the measurable interfer- 
ometer output. The probability P\ of finding the atom 
in the upper energy eigenstate at any time is given by 

Pi = \{M^)\\ (20) 

and this will remain zero unless an asymmetry is present 
together with suitably short splitting and recombining 
stages for the interferometer process - so that transitions 
occur between \cj) } and \c/>i) due to the presence of V as - 
We find that 

Pi = l(l + ci z V-a x Vl-V 2 ) 



At the beginning and end of the interferometer pro- 
cess V <C 1 and hence the probability Pi only depends 
on the x component of the Bloch vector. The probability 
Pi (T) thus depends on how this component has changed 
from its initial value of 1. We can therefore describe the 
dynamical behavior of the single atom interferometer in 
terms of the evolution of the Bloch vector during the 
splitting, holding and recombining stages. At the start 
of the process the Bloch and torque vectors are antiparal- 
lcl [Fig.Qfa)] and approximately aligned with the x axis. 
For small values of the splitting parameter (3 the abso- 
lute value of the torque vector is mainly determined by 
the symmetric energy gap Ao (Fig. ^) and its direction 
remains along the -x direction. The position of the Bloch 
vector remains mostly along the +x direction [Fig. ffib)] 
during early stages of the splitting process. When the 
splitting parameter is increased further the decreasing 
energy gap Ao becomes comparable with and later much 
smaller than the asymmetry quantity V as . As a result 
the torque vector rotates in a x-z plane until it is aligned 
along the z direction [Q « (0, 0, V^ s )]. It is important 
to make this change non-adiabatically so that the Bloch 
vector does not follow the torque rotation. If the Bloch 
vector were to follow the changes of the torque vector adi- 
abatically the atom will always stay in the ground state 
and no interference would be observed. 

During the holding stage the torque vector is con- 
stant and the Bloch vector precesses around the torque 
vector with a constant angular velocity V as , and hence 
both the x and y components oscillate with a period 
2ir/Vas [Fig- Etc)]. In an ideal double- well interferom- 
eter the splitting and recombination stages are short and 
the value of the x component does not change much dur- 
ing these stages, so that a x (T) (which defines the final 
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excited state population) is basically given by its value 
at the end of the holding period. The simple behav- 
ior during the holding period indicates that the excited 
state population would have a period 2-k /V as considered 
as a function of holding time. A similar description in 
terms of the evolution of a Bloch vector also applies to 
the scheme described in Ref. , though the dynamical 
behavior of the Bloch vector is different. 

The behavior of the interferometer may also be de- 
scribed in terms of time-dependent states \L), \R), which 
during the holding period represent atoms localized in 
the left and right wells. The interferometer process in- 
volves the transition |5(0)) — > \AS(T)), which involves 
two pathways \S(0)) — ► \L(T/2)) — ► \AS(T)) and 
|S(0)) — ► \R{T/2)) — > \AS(T)), involving two possi- 
ble localized intermediate states associated with the left 
or right wells. The overall transition amplitude is the 
sum of amplitudes for the two pathways, and depending 
on the relative phase between these amplitudes cither 
constructive or destructive interference may occur. For 
maximum contrast it is desirable that the magnitudes of 
the two partial amplitudes be equal, so that during the 
holding period the populations of the left and right well 
states should be about the same. After optimal splitting 
the z component of the Bloch vector a z should be kept 
close to zero during the holding period, however a phase 
difference between the localized states accumulates. Only 
at the end of the recombination stage this phase is trans- 
lated into the population of the excited state. 



V. RESULTS OF NUMERICAL SIMULATIONS 

We studied the evolution of a Bloch vector and the 
population of the excited state by solving Equations l|17|) 
numerically. The splitting parameter (3 is changed lin- 
early from zero up to a maximum /3 max during the split- 
ting period. It is then held constant at /3 max during the 
holding period, and then changed linearly to zero during 
the recombination period. The dynamical behavior of 
the Bloch vector components is shown in Fig. 0Ja) along 
with the time dependence of the asymmetry parameter 
V = V as /A, the splitting parameter j3 and the excited 
state population Pi [Fig. 0(b)] . The parameters used 
are vZ = 0.02 x , /3 max = 12.5 and T s = 20, T h = 20, 
T r = 20 in dimcnsionless harmonic oscillator units. Here 
we observe complex oscillatory behavior for the x and y 
components of the Bloch vector which occurs during the 
splitting and merging stages. During the holding stage 
they exhibit simple periodic variations with frequency 
V as = 0.2. At the same time the z component develops 
a small negative value during splitting and increases the 
absolute value even further during merging. The x com- 
ponent reaches a negative value of -0.9 at the end of the 
process. This corresponds to an excited state population 
of 0.95 and represents a case of constructive interference. 

By monitoring the behavior of P\ during the interfer- 
ometric process we can see when non-adiabatic evolu- 




time 



FIG. 4: (a) Time evolution of the Bloch vector components 
(j x (solid line), a y (dashed line) and a z (dotted line) for 
T s = T h = T r = 20 and /3 max = 12.5; (b) Time evolution of 
the first excited state population Pi (dotted line), the asym- 
metry parameter V (dashed line) and the splitting parameter 
P/Pmax (solid line). 

tion occurs. The population Pi changes from to 0.47 
[Fig. IHb)] at the beginning of the splitting process and 
does not reach the optimal value 1/2 as a result of the 
non-zero z component of the Bloch vector. The variable 
Pi does not change during the adiabatic precession of the 
Bloch vector around the torque vector during the rest of 
the splitting, holding and the beginning of recombining 
stages. It again exhibits drastic changes in a short pe- 
riod during the re-merging when the torque vector rotates 
rapidly and the Larmor frequency is relatively small. 

It is tempting to limit the evolution of the Bloch vec- 
tor and the relevant phase accumulation during the split- 
ting and merging stages by making these stages shorter. 
However this can lead to excitations of higher excited 
states. We have compared outcomes of the Bloch vec- 
tor model with the results of the numerical solution of a 
multi-state Schrodinger equation using the XMDS code 
[30| . The behavior of the excited state population Pi (T) 
at the end of the interferometer process as a function of 
the holding period Th is shown in Fig. El for the parame- 
ters KT = 0.02 x , /? max = 12.5, T s = T r = 5 [Fig. Eta)], 
T s = T r = 20 [Fig. 13b)] and T s = T r = 200 [Fig. Etc)]. 
In all cases the sinusoidal behavior of the excited state 
population as a function of the holding period can be 
seen. Situations ranging from complete destructive in- 
terference to perfect constructive interference are both 
present. For short splitting times [Fig. Eta)] we observed 
a discrepancy between the two models. Multi-state nu- 
merical simulations indicate the presence of populated 
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FIG. 5: Dependence of the first excited state population Pi 
at the end of the interferometric process on the duration of 
the holding stage Th for various durations of the splitting 
and recombining stages T s — T r = 5 (a), 20 (b) and 200 (c). 
Results of the Bloch vector model are represented by solid 
lines and outcomes of full numerical simulations are presented 
by circles. 



higher energy states which the Bloch vector model ig- 
nores. The full numerical calculations show an irregular 
high frequency modulation of the fundamental frequency 
signal and a reduced maximum population of the first 
excited state. 

For splitting and merging times T s = T r = 20 
[Fig-OJk)] both models show excellent agreement indicat- 
ing a simple sinusoidal variation of the first excited state 
population with holding time. This simple behavior is 
also observed for long splitting time [Fig. EJc)] but with 
significantly reduced amplitude of the oscillations. The 
reduced fringing is attributed to the onset of adiabatic 
following of the Bloch vector during splitting and recom- 
bination which is shown by the presence of a non-zero 
a z component [Fig.^Ja)]. We noticed that our numerical 
solution of the Bloch equation is robust with regard to 
the variations of the signal but is fragile regarding the 
phase. The error was accumulated during the splitting 
stage as a result of V as ^ in a merged trap and will 
scale with the splitting time. 

In the asymmetric double-well potential the ground 
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FIG. 6: Dependence of the filling factor F on the duration of 
the splitting stage T 3 for /3 max = 12.5 and various values of 
asymmetry VT S = 0.01 x (a), 0.02 x (b), 0.05 x (c) and 0.1 x(d). 

state eigenfunction will predominantly occupy the lower 
well, and the excited state eigenfunction will be local- 
ized in the upper well [Fig. |2{c)]. In the slow splitting 
regime the onset of the adiabatic evolution will lead to 
the unbalanced distribution of the atomic wavefunction 
between the wells, which in turn leads to a reduction in 
the measured signal. In application to interferometry it 
can be seen as intrinsic which-way information when the 
atom will predominantly follow one path after the split- 
ting. In general 

|*) = a|to> + &|&>+|&>, (22) 

where | 4>i ) is a linear combination of all other excited 
states. We define a filling factor 

F = 2ab, (23) 

which will describe the balance of ground and excited 
states populations. The dependence of the filling fac- 
tor on splitting time is shown in Fig. El for a splitting 
of (3 = 12.5 and different asymmetries. The results of 
the Bloch model (dotted line) and the multi-state nu- 
merical simulations (solid line) show good agreement for 
splitting times T s > 20. For shorter splitting stages the 
two mode approximation fails and excitations into higher 
modes take place. In the case of the high asymmetry 
V as = 0.1 x [Fig. Eld)] we observe a significant devia- 
tion between outcomes of two models. For large values 
of the asymmetry frequency V as it is impossible to adi- 
abatically isolate two lower states and the transitions to 
higher states have to be taken into account. 

VI. DISCUSSION AND CONCLUSIONS 

We have applied a Bloch vector model to describe the 
quantum state interference of a single-atom wavefunc- 
tion in a time- variable asymmetric double-well potential. 
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The probability of finding the atom in the first excited 
state is closely associated with the magnitude of the spa- 
tially dependent external potential and could be used as a 
measure of the applied asymmetry. Transitions between 
ground and excited states occur during the splitting and 
recombination stages. Larmor precession of the Bloch 
vector during the holding stage is induced by the asym- 
metry, will effect an interferometric phase and determine 
the final value of the excited state population. The evolu- 
tion of the Bloch vector during the splitting and merging 
stages is also important because it will affect the mea- 
surable probability Pi. We have also shown that spe- 
cial requirements apply to the duration of the splitting 
and merging stages in order to avoid excitation of higher 
modes for short times and partial adiabatic following if 
the splitting is too long. Both these effects lead to a de- 
crease of the measured signal. Interestingly enough they 
do not affect the contrast of the interference fringes if the 
first excited state is not initially populated. 

Adiabatic evolution of the Bloch vector can offer a new 
way to measure the first excited state population after 
the double-well interferometric process. We have already 
mentioned that in the far-split regime the excited state 
wavefunction does not overlap with the ground state 
wavefunction and will predominantly occupy the higher 
energy well [Fig.^c)]. If at the end of the non-adiabatic 
splitting, phase evolution and non-adiabatic recombina- 
tion process we also add an additional stage of adiabatic 
splitting in a known asymmetrical potential, then the 
wavefunctions of the two states will be spatially sepa- 
rated. For recording the output Pi we now simply mea- 
sure the population of the higher energy well. To shorten 
the adiabatic evolution time we need to apply the highest 
available asymmetry (Fig.[SJ). 
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APPENDIX: DERIVATION OF BLOCH 
EQUATIONS 



U (t)\^>(0)). Operators in the Hciscnbcrg and Schrodinger 
pictures are related via U as = (Oy {h) s (U) . The 

expectation values of operators in the two pictures are 
related as (fi) = (*(i)|(fi) s |*(i)) = (tf (0)|(fi) ff |*(0)). 

The equation of motion for the Bloch vector compo- 
nents can be derived using the Heisenberg picture via 



(A.l) 



= wo)i|(^n*(o)>, 

where the Heisenberg equation of motion for the Pauli 
spin operator in dimensionless units is 



H 



The derivation involves the use of the commutation 
rules for the Pauli spin operators. For the first term, we 
have after substituting for H from Eq. I|14(l 



—i 






MM 


T 






(n 


x0) 



b=x,y,z 



(A. 



Hence the contribution from the first term in Eq. (| A. 2|) 
is given by 



For the contribution from the second term in Eq. 
(|A.2|I . we may first write al = k ab \ A ) ( b V 

A,B=L,R 

where the K\ B are time independent coefficients that 
can be read from Eqs. and then 



a 

dl {<7a) 



H 



E K AB 

A,B=L,R 



(||A»<B| + |A)(|<B|) 



H 



(A.5) 



The state vectors at time t and at time are re- 
lated via the unitary evolution operator U(t) as = 



Using Eq. (|12|l for the state vector and reverting to 
the Schrodinger picture we have 



J 



(*(0)\(d t (^ a )f\*(0)h = (*(t)\i £ K a AB [(d t \A))(B\ + \A)(d t (B\)}} |*(./l> 

I A,B=L,R 

= E K% B [C* D C B (D\(d t \A)) + C* A C D (d t (B\)\D)} 

A,B,D=L,R 



(A.6) 
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To evaluate this result, a consideration of the four 
quantities (A\ (d t \B)), where (A, B = L, R) is required, 
noting also that {d t {B\)\A) = {{A\ {d t \B)))* . These 
four quantities can be expressed in terms of related ma- 
trix elements in the symmetric basis (Si \ (dt\Sj)), where 
(i, j = 0, 1). Note |Sb) = \S) and \Si) = \AS). 

Using the normalisation and reality property, we first 
show that the diagonal terms (Si \ (dt\Si)) are zero. For 
the off-diagonal terms (d t \Sj)), these are zero be- 
cause \Si) and d t \Sj) have opposite symmetry. From 
these considerations we find that all matrix elements 
(A\ (dt\B)), where (A, B — L, R), are zero. Hence the 



second contribution to the equation of motion for the 
Bloch vector component is zero 

<*(0)||(^)*|*(0)> 2 = 0. (A.7) 
Combining both contributions we find that 

^ a =^,a) a . (A.8) 

Thus the Bloch equations can be expressed in vector form 
as in Ea. H18fl and in detail as in Eqs. (|17fl . 
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